Anisotropic Hydrodynamics, Bulk Viscosities and R-Modes of Strange Quark Stars 

with Strong Magnetic Fields 



Xu-Guang Huang 1,2,3,4 , Mei Huang 3 , Dirk H. Rischke 1,2 , and Armen Sedrakian 2 
1 Frankfurt Institute for Advanced Studies, D-60438 Frankfurt am Main, Germany 
Institut fur Theoretische Physik, J. W. Goethe-Universitdt, D-60438 Frankfurt am Main, Germany 

China 



o 

(N 
(N 

9 L1, 
6 



> 
m 
m 



On 
O 



13 



Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100039, 
4 Physics Department, Tsinghua University, Beijing 100084, China 
(Dated: February 26, 2010) 

In strong magnetic fields the transport coefficients of strange quark matter become anisotropic. 
We determine the general form of the complete set of transport coefficients in the presence of a strong 
magnetic field. By using a local linear response method, we calculate explicitly the bulk viscosities 
and Gi transverse and parallel to the B-field respectively, which arise due to the non-leptonic 
weak processes u + s <B> u + d. We find that for magnetic fields B < 10 17 G, the dependence of £j_ 
and £u on the field is weak, and they can be approximated by the bulk viscosity for zero magnetic 
field. For fields B > 10 18 G, the dependence of both £j_ and £y on the field is strong, and they 
exhibit de Haas- van Alphen-type oscillations. With increasing magnetic field, the amplitude of these 
oscillations increases, which eventually leads to negative £j_ in some regions of parameter space. We 
show that the change of sign of C,±_ signals a hydrodynamic instability. As an application, we discuss 
the effects of the new bulk viscosities on the r-mode instability in rotating strange quark stars. We 
find that the instability region in strange quark stars is affected when the magnetic fields exceeds 
the value B = 10 17 G. For fields which are larger by an order of magnitude, the instability region is 
significantly enlarged, making magnetized strange stars more susceptible to r-mode instability than 
their unmagnetized counterparts. 

PACS numbers: 12.38.-t, ll.10.Wx, 12.38.Mh, 21.65. Qr 



I. INTRODUCTION 

Neutron stars provide a natural laboratory to study 
extremely dense matter. In the interiors of such stars, 
the density can reach up to several times the nuclear 
saturation density, n ~ 0.16fm~ 3 . At such high den- 
sities quarks could be squeezed out of nucleons to form 
quark matter [1-3]. The true ground state of dense 
quark matter at high densities and low temperatures 
remains an open problem due to the difficulty of solving 
non-pcrturbativc quantum chromodynamics (QCD). It 
has been suggested that strange quark matter that con- 
sists of comparable numbers of u, d, and s quarks may 
be the stable ground state of normal quark matter [4] . 
This led to the conjecture that the family of compact 
stars may have members consisting entirely of quark 
matter (so-called strange stars) and/or members fea- 
turing quark cores surrounded by a hadronic shell (hy- 
brid stars) [5]. 

Observationally, it is very challenging to distinguish 
the various types of compact objects, such as the 
strange stars, hybrid stars, and ordinary neutron stars. 
Their early cooling behavior is dominated by neutrino 
emission which is a useful probe of the internal compo- 
sition of compact stars. Thus, cooling simulations pro- 
vide an effective test of the nature of compact stars [6- 
16]. However, many theoretical uncertainties and the 
current amount of data on the surface temperatures of 
neutron stars leave sufficient room for speculations [17- 
19]. Another useful avenue for testing the internal 
structure and composition of compact stars is astro- 
seismology, i.e., the study of the phenomena related 
to stellar vibrations [20-25]. In particular, there are 
a number of instabilities which are associated with the 
oscillations of rotating stars. Here we will be concerned 
with the so-called r-mode instability [see Refs. [22, 23] 



for reviews] . This instability is known to limit the an- 
gular velocity of rapidly rotating compact stars. The r- 
mode and related instabilities in rotating neutron stars 
are damped by the shear and bulk viscosities of mat- 
ter, therefore these are important ingredients of theo- 
retically modelling rapidly rotating stars. Such models 
and their microscopic input can then be constrained via 
the observations of rapidly rotating pulsars, such as the 
Crab pulsar and the millisecond pulsars. 

For quark matter in chemical equilibrium, the shear 
viscosity is dominated by strong interactions between 
quarks. The bulk viscosity, however, is dominated 
by flavor-changing weak processes, whereas strong in- 
teractions play a secondary role. For normal (non- 
superconducting) strange quark matter, the bulk vis- 
cosity is dominated by the non-leptonic process [26-30] 

u + s — > u + d, (la) 

u + d — > u + s, (lb) 

since the contributions of the lcptonic processes u + e O 
d + v and u + e <H> ,s + v arc suppressed due to much 
smaller phase spaces. The bulk viscosity of various 
phases of quark matter has been studied extensively, 
see Refs. [24, 26-41]. 

Compact stars are strongly magnetized. Neutron 
star observations indicate that the magnetic field is of 
the order of B ~ 10 12 — 10 13 G at the surface of ordi- 
nary pulsars. Magnetars - strongly magnetized neutron 
stars - may feature even stronger magnetic fields of the 
order of 10 15 - 10 16 G [42-48]. An upper limit on the 
magnetic field can be set through the virial theorem. 
Gravitational equilibrium of stars is compatible with 
magnetic fields of the order of 10 18 - 10 20 G [49-51]. 
In such a strong magnetic field, not only the thermo- 
dynamical but also the hydrodynamical properties of 
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stellar matter will be significantly affected. In partic- 
ular, due to the large magnetization of strange quark 
matter the fluid will be strongly anisotropic in a strong 
magnetic field (we note here that the magnetization 
of ordinary neutron matter is small [52]). Therefore, 
there is need to develop an anisotropic hydrodynamic 
theory to describe strongly magnetized matter in com- 
pact stars. As wc show below, the matter is completely 
described in terms of eight viscosity coefficients, which 
include six shear viscosities and two bulk viscosities. 

In this paper, we will carry out a theoretical study of 
the anisotropic hydrodynamics of magnetized strange 
quark matter and will calculate the two bulk viscosities. 
We will also discuss the implications of the anisotropic 
bulk viscosities on the r-mode instability in rotating 
quark stars. 

The paper is organized as follows. The formalism 
of anisotropic hydrodynamics for magnetized strange 
quark matter is developed in Sec. II. In Sec. Ill we 
apply the local linear response method to derive explicit 
expressions for bulk viscosities. The stability of the 
fluid under strong magnetic field is analyzed in Sec. IV. 
Section V contains our numerical results for the bulk 
viscosities. The damping of the r-modc instability in 
rotating quark stars by the bulk viscosity is studied in 
Sec. VI. Section VII contains our summary. We use 
natural units H = ks = c = 1. The metric tensor is 
g 1 *" = diag(l, -I, -I, -I). We will use the SI system 
of units in our equations involving clcctromagnctism, 
however we will quote the strength of the magnetic field 
in CGS units (Gauss), as is common in the literature 
on compact stars. 

II. ANISOTROPIC HYDRODYNAMICS 

A. Ideal Hydrodynamics 

Hydrodynamics arises as an effective theory valid 
in the long-wavelength, low-frequency limit where the 
energy- momentum tensor T"", the conserved baryon 
current rig, the conserved electric current n£, and the 
entropy density flux s^, etc., are expanded in terms of 
gradients of the 4-velocity and the thermodynamic 
parameters of the system, such as the temperature T, 
baryon chemical potential fj. b , etc.. The hydrodynamic 
equations can be expressed as conservation laws for the 
total energy-momentum tensor T^ u , as well as baryon 
and electric currents, n B and n%. The zeroth order 
terms in the expansion correspond to an ideal fluid and 
we shall use the index to label them. In the presence 
of an electromagnetic field, the zeroth-order terms can 
be generally written as [53, 54], 

rpflU rpflU , rpflU 

= eu»u v - PA"" -^(M" A f A " + M vX F x ") , 

s% = su», (2) 

where e, P, ns, n e , and s are the local energy den- 
sity, thermodynamic pressure, baryon number density, 



electric charge density, and entropy density, respec- 
tively measured in the rest frame of the fluid. A"" = 
' — u^u v is the projector on the directions orthogonal 
to vt 1 . 

Here = -F» X F\+ g^FP a F pa /4: is the energy- 
momentum tensor of the electromagnetic field. F^" is 
the field-strength tensor which can be decomposed into 
components parallel and perpendicular to as 

F" v = F^uxv? - F vX u x u tl + A^P^A/ 

= E»u v - E v u^ + l e ^ al3 (u a B p - u l3 B a ) , 

(3) 

where in the second line we have introduced the 4- 
vectors P" = F^ v u v and = ^ vaf} F va u p /2 with 
e fivaP being the totally anti-symmetric Levi-Civita ten- 
sor. In the rest frame of the fluid, = (I, 0), we have 
E° = B° = 0, E l = F M and B l = ~e ljk F jk /2, which 
are precisely the electric and magnetic fields in this 
frame. Therefore, P" and P M are nothing but the elec- 
tric and magnetic fields measured in the frame where 
the fluid moves with a velocity w M . 

The antisymmetric tensor M^ v is the polarization 
tensor which describes the response to the applied field 
strength F^ v . For example, if f2 is the thermodynamic 
potential of the system, M^ LV = — dVl/dF pu . For later 
use, we also define the in-medium field strength tensor 
H nv _ F ^ _ M p,v_ In ana i g y to F**" we can decom- 
pose M' 1 " and H» v as 

M» v = (PV - P^u u ) + I e ^ Q ' 3 {Mpua - M a up) , 

(4) 

with P' 1 = -M^Uv, M" = e^ afs M ua u p /2, = 
H^u u and = e^"- 9 H ua u fi /2. 

In the rest frame of the fluid, the non-trivial 
components of these tensors are (F 10 , F 20 , F 30 ) = 
E, (P 32 , P 13 , F 21 ) = B, (M 10 ,M 20 ,M 30 ) = P. 
(M 32 ,M 13 ,M 21 ) = M, (P 10 ,P 20 ,P 30 ) = D, and 
(P 32 , H 13 , H 21 ) = H. Here P and M are the electric 
polarization vector and magnetization vector, respec- 
tively. In the linear approximation they are related 
to the fields E and B by P = x e E and M = x m B, 
with Xe and Xm being the electric and magnetic sus- 
ceptibilities. The 4- vectors P", P M , • ■ ■ are all space- 
like, P M u M = 0, B^u^ — 0, ■ • • , and normalized as 
E^Ef, = -E 2 1 B^B p = -P 2 , • • • , where E = |E| and 
B = |B|. 

Since the electric field is much weaker than the mag- 
netic field in the interior of a neutron star, we will ne- 
glect it in most of the following discussion. Upon intro- 
ducing the 4- vector 6" = P M /P, which is parallel to P M 
and is normalized by the condition b^b^ = — 1, and the 
antisymmetric tensor b^ v = e tiVa ^b a up 1 we can write 

= -Mb"", 

= -Hb^, (5) 
with M = |M| and H = |H|. 
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The Maxwell equation e^ al3 dpF v 
form 



takes the 



V x (v x B), 



8 V (B»u' y - B»u») = 0. (6) 

Its non-relativistic form, which is known as the induc- 
tion equation, is given by 

dB 

~dt 

VB = 0, (7) 

where v is the 3-velocity of the fluid. Contracting 
Eq. (6) with gives 

e + D^B-vTWd^bv =0, (8) 



where 9 



d^u^ and D 



u^d u . The second Maxwell 



equation can be written as 



(9) 



V x (H - D x v) 



H x d t v 



(10) 



whose non-relativistic form is 

V D H V xv 
<9D 

~dt 

where n° is the electric charge density and n e is the cor- 
responding current. It is useful to re- write the energy- 
momentum tensor in the following form [55-57], 

>i*« v - PiS"" + P\\b^b v , 



J F0 — 



£vru 
1 



?Zl = ^(uW-E^-bn"), 



(11) 

where "E^" is the projection tensor on the direction per- 
pendicular to both and 6 M . We have defined the 
transverse and longitudinal pressures Pj_ = P — MB 
and P|| = P relative to 6 M ; here P is the thermody- 
namic pressure. In the absence of a magnetic field, the 
fluid is isotropic and P± = P\\ = P. In the local rest 
frame of fluid, we have b^ = (0, 0, 0, 1) (without loss of 
generality, we choose the z-axis along the direction of 
the magnetic field), hence the electromagnetic tensor 
takes the usual form, while Tpp = diag(e, P±, P±,P\\). 

Next we would like to check the consistency of the 
terms that appear in Tpp with the formulae of standard 
thermodynamics involving electromagnetic fields. By 
using the thermodynamic relation 



e = Ts + [IbUb + H e n e - 
and the conservation equations for n' B0 , it e0 , 



(12) 



< , and s£ 

in ideal hydrodynamics, one can show that the hy- 
drodynamic equation u^d^T^ = together with the 
Maxwell equation (8) implies 



De = TDs + n B Dn B + /J, e Dn e - AIDB, 



(13) 



which is consistent with the standard thermodynamic 
relation 

de = Tds + [iQdnB + /J, e dn e — MdB. (14) 

One should note that the potential energy —MB has al- 
ready been included in our definition of e. Otherwise, 
new terms -MB, -D(MB) and —d{MB) should be 
added to the left-hand sides of Eq. (12), Eq. (13), and 
Eq. (14), respectively. Thus, we conclude that our hy- 
drodynamical equations are consistent with well-known 
thermodynamic relations. 



B. Navier-Stokes-Fourier-Ohm Theory 

By keeping the first-order terms of the derivative ex- 
pansion of conserved quantities one obtains the Navier- 
Stokcs-Fourier-Ohm theory. In this theory, T^ v . n g , 
n%, and can be generally expressed as 



T£ v + h"u v + h v u >1 + T ftv , 



115) 



where ,t^ v ,j B ,j^ , and jff are the dissipative fluxes. 
They all arc orthogonal to u M ; this reflects the fact 
that the dissipation in the fluid should be spatial. We 
shall assume that j£ can be expressed as a linear com- 
bination of h^jjg, and j£ [58, 59]. This allows us to 
incorporate the fact that the entropy flux is determined 
by the energy-momentum and baryon number diffusion 
fluxes. Thus, 



7^ 



ubJb 



(16) 



with the coefficients 7, a ei and as being functions of 
thermodynamic variables. 

Next, the hydrodynamic equations are specified by 
utilizing the conservation laws of the total energy- 
momentum , the baryon number density flow n B , 
electric current n£, and the second law of thermody- 
namics, 



d^ v = 0, 



- 0, 



V = 

Td^ > 0. 



(17) 



To discuss the dissipative parts, let us first define 
the 4- velocity u^, since it is not unique when energy 
exchange by thermal conduction is allowed for. We will 
use the Landau-Lifshitz frame in which is chosen to 
be parallel to the energy density flow, so that = 0. 
Upon projecting the first equation of Eq. (17) on u v 
and after some straightforward manipulations, we find 

(e + P)0 + Ds- r^dpUv + MDB = j*u"F vX . (18) 

Combining Eq. (18), Eq. (12), and the second equation 
in Eq. (17), we arrive at 



TV 



+( Me - Ta e ) V - XiTVpOe + (19) 

where V,, = A Mv d" and w> w = \ (W + Vu* 4 ). For 
a thermodynamically and hydrodynamically stable sys- 
tem, Eq. (19) should be non-negative. This implies 



W a fty 



' U a e + E u ), (20) 

where f3 = 1/T, rj^ap is the rank- four tensor of vis- 
cosity coefficients, and and arc thermal and 



OLB 


= PftB, 


a e 


= Pfl>e, 




— yMvafi 




= 


ft 
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electrical conductivity tensors with respect to the diffu- 
sion fluxes of baryon number density and electric charge 
density. By definition, r]^ va ^ is symmetric in the pairs 
of indices a, (3 and fj., v. It necessarily satisfies the con- 
dition r]^ a f s {B a ) = ifP^(-B a ), which is Onsager's 
symmetry principle for transport coefficients. Similarly, 
the tensors k' 1 " and should satisfy the conditions 



(— j3 ) and a^{B x ) = cr^(-B A ). Fur- 



thermore, all the tensors of transport coefficients rfvafi ; 
k^ v , and cr^ v must be orthogonal to by definition. 

As we have seen, the appearance of the magnetic 
field makes the system anisotropic. Such anisotropy 
is specified by the vector 6 M , so that the tensors 7/ ll/Q ^, 
a^ v , and k^ v should be in general expressed in terms of 
u M , 6 M , g^ v ., and V*" . All independent irreducible tensor 
combinations having the symmetry of r]^ 01 ^ and which 
are orthogonal to u M are [60] (see also Appendix A) 



i) A^A Q/3 , 

ii) A^A^ + A^A" a , 

iii) A^b a b f3 + A af3 b^b u , 

iv) lfb v b a ^, 

v) A>* a b u b f3 + A^<V& Q + A^b v b a + A ua b^b f3 , 

vi) A^ a b" p + A v ^b l " x + A^b ua + A" a b^, 

vii) b^ a b v b' 3 + b^b^b" + b^b v b a + b va b»b' 3 , 



b^b 



a-LvB 



b^b va . 



(21) 



All independent irreducible tensor combinations having 
the symmetry of and a^ v and which are orthogonal 
to are 



(i) A"", 

(ii) V*V, 

(iii) b" v . 



(22) 



In accordance with the number of tensors (21) and (22), 
a fluid in a magnetic field in general has eight indepen- 
dent viscosity coefficients, three independent thermal 
conduction coefficients and three independent electrical 
conductivities. They may be defined as the coefficients 
in the following decompositions for the viscous stress 
tensor, heat flux, and electric charge flux 



t" u = 2 m {w^ - A^e/3) + 771 (A"" - l E J i ~ \ ( n - 2??2 ( b ^ abP + b u ^ a bP) w a p 
- m {2b^b v < 3 w afi - E^wJ ~ ~» a Wo ?) - 2,74 {^ a b uP + Z» a b^) w a p 

+ 2 V5 (b^Vti 3 + b va b^) w ali + ^(±Z^(f> + H\\VV(p, (23) 

j% = nTVa B - mWVTVvCXB - K 2 b» v TV v a B , (24) 

ft = a(\7^a e + E^)-a 1 b^b u (V f "a e + E^)~<7 2 b fny (\7^a e + E f "), (25) 

I 



where <fi = ^" v w pv , f> = b t± b l, w fJ , ll and r* 1 " is con- 
structed so that the 77's are the coefficients of its trace- 
less parts, i.e., they can be regarded as shear viscosi- 
ties; C's arc the coefficients of the parts with non-zero 
trace and can be considered as bulk viscosities. The 



k's and er's are thermal and electrical conductivities, 
respectively. 

Now the divergence of entropy density flux (19) can 
be explicitly written as 



Td^ = 2770 Lr - ^a^ («v - l A ^ e ) +m(e- §0) 

+ 27 72 (Wbpvf" - b v b pW P») (b^b p w pu - b u b p w pp ) 
+ V3 {b^w" - b v "w/) - b up w^) 

- a(TWa e + E»)(TV > e + E p ) + a^TWV > e + E^brf. (26) 

I 



One should note that the terms corresponding to the transport coefficients 7/4,775, k 2 , and a 2 in Eqs. (23)- 
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(25) do not contribute to the divergence of the entropy 
density flux. For stable systems, all the other trans- 
port coefficients must be positive definite according to 
the second law of thermodynamics. In Sec. IV we will 
demonstrate explicitly that negative bulk viscosities C_l 



or/and £y indeed cause an instability in the hydrody- 
namic evolution of strange stars. 

To conclude this section, we compare our definition of 
the viscosity coefficients in Eq. (23) with the definition 
given in Ref. [60] for non-relativistic fluid, which reads 



Tij = 2f) (wij - SijO/S) + C,8ij9 

+ f)x{2wij - 5^6 + SijWkibkbi - 2w ik bkbj - 2w jk b k bi + bibjd + bibjWkibkbi) 

+ 2fj 2 (w lk b k b : j + Wjkbkbi - 2b l b J w k ib k bi) + fj 3 (w lk b jk + w jk b lk ~ w k ib lk bjbi - w k ib jk bibi) 

+ 2fj4(wkibubjbk+Wktbjibibk)+(i(SijWktbkbi+bibj9), (27) 

where fey = £ij k b k and the remaining notations are self-explanatory. 

Our viscosity coefficients in Eq. (23) are related to the coefficients in Eq. (27) by 

vo = fj + fji, in = | (fji + |Ci - |c) , m = m - fji, ??4 = 5773, 

vn=m, C± = C + iCi. Cn=C + |Ci- (28) 

In Ref. [60] there is no term that corresponds to our rf 3 . The reason is that Ref. [60] considers the combination of 
vectors (viii) in Eq. (21) as dependent on the others; in the Appendix we will show that (at least for relativistic 
fluids) all the combinations (i)-(viii) are linearly independent. Note that the transport coefficients in Eq. (26) 
appear as prefactors of quadratic forms, therefore the second law of thermodynamics requires that these coefficients 
must be positive definite for stable ensembles. This is not manifest in Eq. (27). 



III. BULK VISCOSITIES 

The typical oscillation frequency of neutron stars is 
of the order of magnitude of the rotation frequency, 
Is -1 < u < 10 3 s _1 . The most important microscopic 
processes which dissipate energy on the corresponding 
timescales are the weak processes. 

The compression and expansion of strange quark 
matter with nonzero strange quark mass will drive the 
system out of equilibrium. The processes (la) and (lb) 
are the most efficient microscopic processes that restore 
local chemical equilibrium. Therefore, the bulk viscosi- 
ties are determined mainly by the processes (la) and 
(lb). In this section we will derive analytical expres- 
sions for the bulk viscosities (± and (\\ [87]. 

Let us imagine an isotropic flow v(t) <~ e 1 "* which 
characterizes the stellar oscillation. If there are no dis- 
sipative processes, such an oscillation will drive the sys- 
tem from one instantaneous equilibrium state to an- 
other instantaneous equilibrium state. The appearance 
of dissipation changes the picture: during the oscil- 
lations the thermodynamic quantities will differ from 
their equilibrium values. Let us explore how the ther- 
modynamic quantities evolve during the flow oscilla- 
tion. 

In general, we can write the change of baryon density 
n B = (n u + n c i + n s )/3 induced by the oscillation of the 
fluid as 



n B {t) 
Sn B 



n B0 - 

r eq 



5n B (t), 
- 6n' B , 



(29) 



where n B o is the static (time-independent) equilibrium 
value, Sn C g denotes the equilibrium value shift from 
n B o due to the volume change and Sn' B denotes the in- 



stantaneous departure from the equilibrium value. Be- 
cause processes (la) and (lb) conserve baryon number, 
8n' g (t) can be set to zero, if we neglect other micro- 
scopic processes. Then Sn B can be determined through 
the continuity equation of ideal hydrodynamics, 



5n B (t) = —t 



(30) 



Since processes (la) and (lb) also conserve the sum 
rid + n s , a similar argument leads to the relation 



Srid + Sn s 



rido + n s0 



0. 



(31) 



When the system is driven out of chemical equilib- 
rium, the chemical potential of the s quark will be 
slightly different from that of the d quark. Let us de- 
note this difference by Sfi = fi s — fid = Sfi s — 5 fid, with 
Sfif being the deviation of fij from its static equilibrium 
value. Up to linear order in the deviation we find 



6fi{t) 



dn s J s \dnd/ 



5n d , 



(32) 



where rif denotes the number density of quarks of fla- 
vor / and the subscript indicates that the quantity in 
the bracket is computed in static equilibrium state. 5n s 
and Srid are the deviations of s-quark and d-quark den- 
sities from their static equilibrium value. In the final 
expressions they should be functions of 0. 

The instantaneous departure from equilibrium is re- 
stored by the weak processes (la) and (lb). Adopting 
the linear approximation, this can be described by 



r s = xSfi, a > 0, 



(33) 
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where Td and T s are the rates of processes (la) and 
(lb), respectively. If the weak processes are turned off, 
one should have 



i c f q = -ri/00 = ri/o 



riBO 



(34) 



where the dot denotes the time derivative. After turn- 
ing on the weak processes, we have 



8n u = 8n B , 

8n B 



Sn d 



5n s 



n d0 - 



riBo 
5n B 



XS^(t), 



n s0 \8fi(t). 



n B Q 



(35) 



This system of coupled linear first-order equations is 
closed by substituting Eqs. (30)-(32). It is then easy to 
obtain the solution, 



Sn u 
Sn d = - 
Sn s = - 



-n u o- 



iojn dQ + A (d[i s /dn s ) (n d0 + n s0 ) 9 
iuj + XA iu ' 

iojn s0 + A (dfj, d /dn d ) (n d0 + n s0 ) 
iu + XA iuj ' 



where the coefficient A is defined by 



A 



dn d 



(36) 



(37) 



The parallel and transverse components of the pres- 
sure Pi i and Pi can be written as 



P,l = P|| cq + <5P||', 
Pi = P± cq + SP ± ', 



(38) 



and 



p eq ^ n O , 



1/ 

where 




dB 



(39) 



Sn'f = Snf — Sn eq . 



(40) 



The small departure of the magnetic field SB can be 
calculated by the variation of Eq. (8). One finds 



ib - -\—e 

6 ioj 

A direct calculation then gives, 
S< = 0, 
Sn d - 



(41) 



5n' = - 



iuj + XA iu> 
AC,, 9 



IU1 



XA 



iu> 



(42) 



where we introduced the coefficient Cii as 
.'d^ d \ ( dfx 



<9n 5 



Now we obtain 



sp/ = XC W C W/i. o 
"A + AA iw' 



(43) 



(44) 



with C± defined as 

C± — C|| — XB, 



X 



fdM\ (9M\ 
\dn d ) \dn s J 



(45) 



Then the deviation of T^ v from its equilibrium value 
can be written as 

ST^ = -ReSP^'E^ + Red P| ( (46) 

For isotropic flows, we have 

= ( 1 _E» u 6-( ll b ti b v e. (47) 

By comparing the above two expressions, we obtain 



and 



Cii 



AC,, 2 



a; 2 + A 2 A 2 ' 

XC±C\\ 
w 2 + A 2 A 2 ' 



(48) 



(49) 



Expressions (48) and (49) show that the bulk viscosities 
(j_ and 0i are functions of the perturbation frequency 
w, the weak rate A and the thermodynamic quantities 
Cii, C±, A. From the derivation above we can convince 
ourselves that these expressions should be valid also 
in the case of color-superconducting matter. For zero 
magnetic field, Eq. (49) reduces to Eq. (48) which, with 
parameters Ch,A, and A taken in the absence of mag- 
netic field gives the expression for the usual bulk vis- 
cosity Co defined in isotropic hydrodynamics. 

Both £_l and 0| attain their maxima in the limit 



SB, of zero frequency, (V 



C||7(aa 2 ), a r 



C||Cj_/(Aj4 2 ); and the maxima are inversely propor- 
tional to the weak interaction rate. At high frequency, 
u ^s> XA, (± and 0| fall off as 1/uj 2 . For practical appli- 
cations to cold strange stars, where the chemical poten- 
tial is much larger than the temperature, the quantities 
C|| , C± , and A can be evaluated in the zero-temperature 
limit. Their dependence on temperature is weak. Con- 
trary to this, the coefficient A depends strongly on tem- 
perature: for normal quark matter, A has a power-law 
dependence on T (see Sec.V); for a fully paired color- 
superconducting phase, the weak rate is exponentially 
suppressed by a Boltzmann factor e~ A / T with A be- 
ing the superconducting gap. Consequently, £j_ and 0| 
depend exponentially on T [24, 32-37, 40, 41]. 

Before we find the numerical values for the bulk vis- 
cosities, we first need to analyze the stability of magne- 
tized strange quark matter. We observe that according 
to Eq. (49) negative values of (± arc a priori not ex- 
cluded. In the following we analyze the consequences 
and implications of negative £_l on the stability of the 
system. 
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IV. STABILITY ANALYSIS 
A. Mechanical Stability 

Stable equilibrium in a self-gravitating fluid, such as 
in strange quark stars, is attained through the balance 
of gravity and pressure. The gravitational equilibrium 
requires that both components of the pressure Pii and 
P± should be positive (otherwise the star will undergo a 
gravitational collapse). At zero temperature, the one- 
loop thermodynamic pressure P — ThiZ/V of non- 
interacting strange quark matter, where Z is the grand 
partition function, is given by 



P 



E 



N c q f B 
4tt 2 



E Vn 



Hfy — to 2 — 2nqjB 



-{m 2 f + 2q f Bn) In 



fi 2 — to 2 — 2qfBn 



m 2 + 2qfBn 



(50) 



where qf is the absolute value of electric charge, fif and 
to/ are the chemical potential and the mass of quark 
of flavor /, n labels the Landau levels, v n = 2 — S 0n 
is the degree of degeneracy of each Landau level and 
n max = IntKA 4 / — TO /)/(2<Z/-B)] is the highest Landau 
level for quarks of flavor /. By differentiating Eq. (50) 
with respect to B one can easily get the magnetization 
as 



7V 



EE- 

/ n=0 



' 4tt 2 



Vfyvj ~ m f~ 2n vB 



-(m 2 f + AnqfB) In 



(if + 4 / fij — mj — 2nqfB 



to 2 + 2nqfB 



(51) 



In Fig. 1 we illustrate the magnetization as a function 
of magnetic field at zero temperature. The parameters 
are chosen as 



Hu = = l-h = 400 MeV, 
m s = 150 MeV, 
m u = rrid = 5 MeV. 



(52) 



On average, the magnetization increases when B grows 
and eventually becomes constant when B > B c , where 



B c ee Max/{(^ - m})/(2q f )}. 



(53) 



However, the detailed structure of the magnetization 
exhibits strong de Haas- van Alphcn oscillations [62]. 
This oscillatory behavior is of the same origin as the 
de Haas-van Alphcn oscillations of the magnetization 
in metals and originates from the quantization of the 
energy levels associated with the orbital motion of 
charged particles in a magnetic field. The irregularity 
of this oscillation shown in Fig. 1 is due to the unequal 
masses and charges of u, d, and s quarks. 

When B > B c , all quarks are confined to their lowest 
Landau level and their transverse motions are frozen. 



In this case, the longitudinal pressure Pu oc B, so the 
magnetization M = dP/dB = P/B is independent of 
B, and the transverse pressure Pj_ = P — MB of the 
system vanishes. This behavior is evident in Fig. 2. 
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FIG. 1: The magnetization of strange quark matter as func- 
tion of the magnetic field B. 
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FIG. 2: (Color online) The parallel P\\ (dashed, black on- 
line) and transverse P± (solid, red online) pressures of 
strange quark matter as functions of the magnetic field B 
in units of the pressure Po for zero magnetic field. 



Thus, we conclude that when all the quarks are con- 
fined to their lowest Landau level, the transverse pres- 
sure vanishes. The system therefore becomes mechan- 
ically unstable and would collapse due to the grav- 
ity [63]. This phenomenon establishes an upper limit 
on the magnetic field sustained by a quark star. Given 
fid ~ 400 MeV, this upper limit is roughly B c ~ 10 20 
G as shown in Fig. 2. 



B. Thermodynamic Stability 

The thermodynamic stability requires that the local 
entropy density should reach its maximum in the cqui- 
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librium state [64] . The total energy density of the fluid 
and the magnetic field is 

B 2 

etotai = Ta + HfUf -P+ — , (54) 

and the corresponding first law of thermodynamics, in 
variational form, is 



fetotai = TSs + HfSrif + HSB, 



(55) 



where H is the strength of the magnetic field and 5 
stands for a small departure of a given quantity from 
its equilibrium value. Varying Eq. (55) on both sides 
and taking into account that etotai, H/, and B arc inde- 
pendent variational variables, one obtains 



5*8 



8n f 5n } 



1 

L 

T 



-5HSB 



8x T x8x, (56) 



where Sx = {5n u , Srid, Sn s , SB) and 



o 





dH 









d^ u 
dB 


d^d 





d Hd 


dn d 


dB 







d[±* 


dn * 

ML 


dB 


dH 


dH 


dn d 


dn s 


dB 



(57) 



Jo 



logB — \x plane in Fig. 3 for our parameters (52). The 
TUW is actually very narrow. One may conjecture that 
such an instability may lead to formation of magnetic 
domains [62]. The presence of possible magnetic do- 
mains in neutron star crusts was discussed in Ref . [65] ; 
furthermore, such a possibility for color-flavor-locked 
quark matter was pointed out in Ref. [66] . We will not 
pursue here the study of domain structure and related 
physics, since among other things, this will require us 
to specify the geometry of the system. 



Mechanically unstable region 




3 18.5 



T=0 
TUW near B" 
TUW near B" 
TUW near B? 
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The thermodynamical stability criteria require that 
5s = 0, 8 2 s < 0, or, equivalently, x is positive defi- 
nite. Taking into account the relation (dH/drif)o = 
— (dfif /dB)o, it is easy to show that these criteria are 
equivalent to the requirement 



/ drif \ ^ 
From Eq. (50) we obtain 



fdM\ 



(58) 



FIG. 3: (Color online) The first three thermodynamically 
unstable windows (TUW) for strange quark matter in the 
log B — Hd plane for T = and fi u = ■ 



C. Hydrodynamic Stability 
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2tt 2 



/ij- — mj — 2nqfB 



;(59) 



then it is evident that the condition (<9n//<9/i/) > 
is always satisfied. However, (dM/dB) is divergent 
when B approaches an^O Landau level for each flavor 
quark from below, 



~dB) 
when 



N c q f B 



(nq f f 



2nq f B 



B^Bl = 



2nqj 



-0 H 



(60) 



This shows that strange quark matter will be thermo- 
dynamically unstable just below each Landau level B^ 
(n ^ 0) for every flavor /. The first three thermo- 
dynamically unstable windows (TUW) associated, re- 
spectively, with Bf,B™, and Bf are illustrated in the 



In this subsection we address the problem of hydro- 
dynamic stability within the theory presented in Sec. II. 
The fluid is said to be stable if it returns to its ini- 
tial state after a transient perturbation. Otherwise, 
i.e., when the perturbation grows and takes the fluid 
into another state, the fluid is unstable. Our particular 
goal here is to determine whether a small, plane-wave 
perturbation around a homogeneous equilibrium state 
grows for nonzero (± and C|| [58]. All the other trans- 
port coefficients are set to zero. For this purpose, it 
is sufficient to solve the hydrodynamic and Maxwell 
equations that are linearized around the homogeneous 
equilibrium state 



d^ST^ = 0, 



d^n% = 0, 
d^n* = 0, 
d^H^ = 0, 
E^dJF af3 = 0, 



(61) 



where 
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= 6T£ + 


0J EM + ( '" i 


OJ F0 


= 5eu>*u" 


-6P ± Z» V + SP^b 


01 EM 


= B5B(u 


V - W v - Wb") + 


St^ 


= UE» V 6 






= Subu^ 


+ nsSu^, 


e 






SH 1 *" 


= SHV™ 


+ HSb^, 


SFV 


= 5Bb^ 






= Sb^b" 4 




5<j> 


= -•sr(a IA 6u v + a v Su ll ), 







Wb v ) + B z (Su^u" + u^Su" - 8Wb v - b»8b u ), 



(62) 

Upon linearizing the normalization conditions u^u^ = = — l^u^b^ = 0, one finds that the perturbed 

variables need to satisfy the constraints 

Su^Uf, = = u^Sbp + Su^bf, = 0. (63) 

In the equations above the perturbations are assumed to have the form SQ = 5Qo cxp(ikx), where SQo is constant, 
and the unperturbed quantities arc independent of space and time. The hydrodynamic and Maxwell equations 
need to be supplemented by an equation of state (EOS) in order to close the system. The linearized EOS is given 
by 

5P = c 2 Se + MSB, (64) 

where c 2 = (dP/de)B is the speed of sound. 

In the most general case Eqs. (61)-(64) constitute 15 independent equations, but the equations associated with 
the conservation of ns and n e are decoupled from the others if the EOS is taken in the form (64). Therefore, we 
are left with 13 equations. We work in the rest frame of the equilibrium fluid, u M = (1, 0, 0, 1) and b^ = (0, 0, 0, 1) 
and choose as independent variables 

8Y, = {6u 1 ,du2,8u 3l 5bo,5b ll Sb 2l de,dP ll ,5P A _ 1 SB,SM,5(j),Sip}. (65) 

(N.B. One can choose other independent variables, but the results do not change). The thirteen linear equations 
can be collected into the following matrix form 

GijSYj = 0. (66) 

The matrix G has the following form 

G = 





fci 















































k 2 


-fcl 























-Hk 


-Hko 








Hk 3 













fci + k 2 


-fcl - fc 2 








h 


h 








—ko 


-ko 























Bki 


Bk 2 























-k 











—hki 


-hk 2 


-hk 3 


-HBk 3 








k 








Bk 











—hko 


—hko 








HBk 3 


HBk 3 








fcl + fc 2 


B(h + fc 2 ) 





-Cx(*i- 


- fca) 








—hko 


— H Bh() 


HBki 


HBk 2 





k 3 





-Bk 3 








C\\h 








1 


i 


















































-1 


1 


M 


B 








iki 


ik 2 





























1 











-ik 3 





























1 




















-c 2 


1 





— M 












(67) 

where h — e + P + HB is the total enthalpy. The exponential plane-wave solutions for frequencies fco and 
wave- vectors k satisfy the dispersion relations given by 



det G = 0. 



(68) 
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For the modes propagating parallel (longitudinal modes) or perpendicular (transverse modes) to the magnetic 
field, Eq. (68) has simple solutions. 

1) Transverse modes, k 3 = 0. There are two types of transverse modes. One is solely determined by the Maxwell 
equations and has the following dispersion relation 

h = ±kj_, (69) 

where k\ = k\ + k\, and describes simply an electromagnetic wave. Another solution has the dispersion relation 

iC ± k 2 ± ± ^(e + P^is + P^kl-kiU 2 



2(e + P||) 



~i— c ^ + wrwr (70) 

where the second approximate relation is valid in the long-wavelength limit. This solution represents a sound 
wave propagating perpendicular to the magnetic field. The speed of this sonic wave is \J{e + P± ) / (e + P\\)c s and 
is smaller than the speed c s of a ordinary sound wave. It is seen that positive £j_ implies dissipation of the sonic 
wave, i.e., a decay of the initial disturbance. We conclude that the fluid flow is stable in the case. However, we 
see that for negative the initial disturbance grows and the fluid is unstable. Thus, we conclude that negative 
transverse bulk viscosity implies hydrodynamic instability via growth of transverse sound waves. 

2) Longitudinal modes, k\ = k2 = 0. We find three types of longitudinal modes. The first one is again the 
electromagnetic wave with the dispersion relation 

k = ±k 3 . (71) 

The second one is a transverse wave oscillating perpendicularly to the magnetic field, but traveling along the 
magnetic field lines. It has the dispersion relation 




fc o = ±W , 5 , D ^ fc 3- (72) 



This mode is the Alfven wave whose speed is equal yJBH/(e + P|| + BH). The third longitudinal mode has the 
following dispersion relation 



K\\kl±\l^e + P\\Yclkl-4^ 



2( £ + Pn) 

^±c s k 3 + ^- y (73) 

This mode represents an ordinary sound wave with dissipation due to the longitudinal bulk viscosity £ii . It is 
obvious that if Q < this mode will not decay, rather grow, thus leading to hydrodynamic instability. 



In the next section, we will show that for certain val- 
ues of the parameters, the transverse bulk viscosity £j_ 
could be indeed negative. We emphasize here that this 
does not imply a violation of the second law of thermo- 
dynamics, rather this manifests a hydrodynamic insta- 
bility of the ground state, i.e., small perturbations will 
take the system via this hydrodynamic instability to a 
new state. A candidate state is the one which has inho- 
mogeneous (domain) structure. Both the structure of 
the new state and the transition from the homogeneous 
to the inhomogeneous state are interesting problems 
which are beyond the scope of this study. However, we 
would like to point out a number analogous cases where 
a negative transport coefficient indicates instability to- 
wards formation of a new state with domain structure. 
One such case is the negative resistivity (also known 
as the Gunn effect) in certain semiconducting mate- 
rials [67, 68]. Another case is the negative (effective) 



shear viscosity, which is extensively studied in the liter- 
ature [69-72] . Finally, negative bulk viscosity has been 
investigated in different contexts in [73, 74]. 



V. RESULTS FOR THE BULK VISCOSITIES 



In order to calculate the bulk viscosities, we need 
to determine the coefficients A, Cii , Cj_ and A. From 
Eqs. (50), (51), and (59) we obtain in a straightforward 
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manner 



/ N c q s B ^ 
1 ~2^~ ^ 



" IT; 
n=o VM: 



- mi 



2nq s B 



+ 



N c q d B 



'A, 



n=0 



Cm = n d0 



— n. 



N c q d B x - 

"2^ V" 



-l 



> — = 

^ VMs - m s - 2nq s B 



Cx = C, - ( dM ; dn <to dM i dn 



i am on d0 _ am on sQ 



B. 



-l 



(74) 



where 



dM 



N, 



If 



'fJ? f -m)- 2nq f B 



n=0 



The rate A of the weak processes (la) and (lb) should 
also be affected by a strong magnetic field. The major 
effect of a magnetic field on A is to modify the phase 
space of weak processes (la) and (lb) [75, 76]. Taking 
this into account, one obtains 



x 



1 




- ™ 2 


- 2nq u B J 




1 


-2nq d B^j 


\/^ 








1 






- TO 2 


- 2nq s B J 



where G 2 = G% sin 2 6 C cos 2 <9 C = 6.46 x 10~ 24 MeV" 4 
is the Fermi constant. 

1) When the magnetic field is much smaller than the 
typical chemical potential, say, q d B <C n d , its effect on 
the bulk viscosities is negligible. For typical parame- 
ters (52), this condition holds up to B ~ 10 17 G. In 
this case, the system is practically isotropic, C|| and C_l 
are effectively degenerate with the isotropic Co > the bulk 
viscosity of unmagnetized matter. The zero-magnetic 
field limit results can be obtained easily by replacing 

X)"=o x QfB ~* %kpf, and B — > 0. The bulk viscosity 
for zero magnetic field Co as function of oscillation fre- 
quency uj for various temperature is shown in Fig. 4. 
The "shoulder" structure and the temperature depen- 
dence of Co are easily understood from Eq. (48) and 
have been widely discussed in the literature [24, 28- 
37, 40, 41]. 

2) When the magnetic field is extremely large, say, 
B ^> B c <~ 10 20 G, for our choice of parameters (52) all 
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FIG. 4: (Color online) The isotropic bulk viscosity Co at 
zero magnetic field as function of the oscillation frequency 
u for fx u = fi d = 400 MeV at T = 0.01 (solid, blue online) 
0.1 (dashed, red online) and 1 (black, dotted online) MeV. 



the quarks are confined in their lowest Landau level. In 
this case we obtain 



A 



2^ 2 



N c q d ^, d B 

»3 



(k Fd + k Fs ), 



/'</ 



c, = 0, 



A = 



and therefore 

a = 

Cn * 



4G 2 g 2 g 2 j? 4 T 

57r 3 ^ 2 fc Fs 



45m 4 ^ 2 fc F 



(77) 



WirG 2 qlB*TZ(k Fs + k Fd y 



(78) 



We used the parameters (52) and assumed physically 
interesting frequencies ui < 10 4 s _1 . The bulk viscosity 
C_l vanishes as a consequence of vanishing Pj_ when 
B > B c . Since C|| is now inversely proportional to B 2 it 
approaches zero for large B. Therefore, both C_i_ and C|| 
are suppressed for large B. In contrary, Py is enhanced 
by the extremely large magnetic field, see Fig. 2. 

3) When the magnetic field is strong, but not strong 
enough to confine all the quarks to their lowest Lan- 
dau level, the situation becomes complicated. For our 
chosen parameters (52), this situation roughly corre- 
sponding to the interval 10 17 G < B < 10 20 G. In this 
finite number of Landau levels is occupied, and 
the essential observation is that Cm and C± can be 
negative. The behaviors of Cm and C± are shown in 
Fig. 5 as functions of B. Let us concentrate on the 
few levels just above the value 10 19 G. When B grows 
passing over B% or B s n for each n. both C\\ and C± 
change their sign. More importantly, they have always 
opposite signs. Therefore, in this region, (± is nega- 
tive which leads to hydrodynamic instability (see the 
analysis in Sec. IV C). 

The numerical values of the bulk viscosities C|| and 
C_l are shown in Fig. 6 as functions of B. The pa- 
rameters are those given in Eq. (52). We also fix the 
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FIG. 5: (Color online) Coefficients C« (dashed, blue on- 
line) and C± (dotted, red online) as functions of B at zero 
temperature. 
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temperature T = 0.1 MeV and oscillation frequency 
uj = 2n x 10 3 s _1 . Both £y and £j_ have "quasi-periodic" 
oscillatory dependence on the magnetic field. The two 
boundaries of each "period" correspond to a pair of 
neighboring B^, f = u,d,s and n = 0, 1,2- ••, and 
hence the period is roughly AB ~ 2qfB 2 /kpj for large 
B. Therefore, on average, the period increases as B 
grows. The amplitude of these oscillations also grows 
with increasing magnetic field until B ~ B c . Thereafter 
all the quarks arc confined to their lowest Landau levels 
and Ci vanishes. From Fig. 6 we see that the magni- 
tudes of C_l and £|| can be 100 to 200 times larger than 
their zero field value £o- Due to the unequal masses and 
charges of u, d and s quarks, (ji and behave very ir- 
regularly. We illustrate the zoomed-in curves around 
B = 10 17 and B = 10 18 G in the sub-panels, which 
look more regular. The quasi-periodic structures are 
more evident in these sub-panels. 

The most unusual feature seen in Fig. 6 is that for 
a wide range of field values, the transverse bulk vis- 
cosity £j_ is negative. Therefore, strange quark matter 
in this region is hydrodynamically unstable. Besides 
this hydrodynamical instability, near each B' nl there is 
a narrow window in which thermodynamical instability 
arises. We depict Q and (± in these unstable regions 
by dashed red curves. The solid blue curves correspond 
to the stable regime. 

The magnetic field in a compact star need not be 
homogeneous and may have a complicated structure 
with poloidal and toroidal components. Furthermore, 
the fields will be functions of position in the star be- 
cause of the density dependence of the parameters of 
the theory. Furthermore, the instabilities, described 
above, may lead to fragmentation of matter and for- 
mation of domain structures, where the regions with 
magnetic fields are separated from those without mag- 
netic field by domain walls. Accordingly, only the aver- 
aged viscosities over some range of magnetic field have 
practical sense for assessing the large-scale behavior of 
matter. Averaging over many oscillation periods in the 
stable region, we find that the averaged values of 0i 
and (j_ are much more regular, with their magnitudes 
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FIG. 6: (Color online) Bulk viscosities Q\ and £j_ scaled by 
the isotropic bulk viscosity £o as functions of the magnetic 
field B at fixed frequency uj — 2tv x 10 3 s -1 and temperature 
T — 0.1 MeV. The dashed red curves correspond to viscosi- 
ties lying in the unstable regions and would be not physi- 
cally reachable. Sub-panels show the amplifications around 
10 17 and 10 18 G. Our parameters are given in Eq. (52). 



restricted from to several (o, see Fig. 7. In obtain- 
ing the curves in Fig. 7, we have eliminated the vis- 
cosities lying in the unstable regime. The solid black 
curves are obtained by averaging over a short period 
A log 10 (S/G) = 0.05. The period was chosen such that 
the most rapid fluctuations are smeared out, but the 
oscillating structures over larger scale are intact. The 
short-dashed red curves correspond to averaging over 
an even longer period, Alog 10 (B/G) = 0.5. The result 
of long-period averaging is that C|| first increases slowly 
and then drops down quickly once B > 10 18 5 G; sim- 
ilarly, (± first slowly decreases and then drops down 
very fast for B ~ 10 18 5 G. Such a dropping behav- 
ior reflects the fact that a large number of quarks are 
beginning to occupy the lowest Landau level. 

Wc note that the appearances of thermodynamic, 
mechanical, and hydrodynamical instabilities arc all in- 
duced by the Landau quantization of the quark levels, 
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FIG. 7: (Color online) The averaged bulk viscosities Q\ and 
£x scaled by the isotropic bulk viscosity £o as functions of 
the magnetic field B at fixed frequency uj — 2ty x 10 3 s _1 
and temperature T = 0.1 MeV. The black-solid curve cor- 
responds to averaging over a short period (Alog 10 (B/G) = 
0.05); while the red short-dashed curve corresponds to av- 
eraging over a long period (Alog 10 (S/G) = 0.5). The vis- 
cosities lying in the unstable regions have been eliminated 
in the averaging. 



i.e., are quantum mechanical in nature. More precisely, 
they are all due to the interplay between the Landau 
levels and the Fermi momentum (reflected in the quan- 
tity B(/). Additionally, the hydrodynamical instabil- 
ity requires that the quark matter is para-magnetized. 
Although we did our analysis by using the free quark 
gas approximation, Eq. (50), it should be valid as long 
as there are sharp Fermi surfaces (low temperature), 
quantized Landau levels (high magnetic field) and para- 
magnetization. The appearances of these instabilities 
are expected to be robust feature for such systems. 

We also checked that if one imposes the neutrality 
condition, i.e., the condition 2n u = rid + n-s, there is 
only minor quantitative change, while the qualitative 
conclusions are almost unchanged. Our choice of chem- 
ical potential p u = pd = 400 MeV roughly corresponds 
to the choice of ~ 4 — 5hq for neutral strange quark 
matter. 



VI. R-MODE INSTABILITY WINDOW 

The purpose of this section is to discuss the damping 
of the r-modes of Newtonian models of strange stars 
by dissipation driven by the bulk viscosities £j_ and £ji . 



As is well known, rotating equilibrium configurations of 
self-gravitating fluids are susceptible to instabilities at 
high rotation rates. Starting from the mass-shedding 
limit and going down with the rotation rate, the first 
instability point corresponds to the dynamical insta- 
bility of the I = 2 and m = 2 mode. This bar-mode 
instability is independent of the dissipative processes 
inside the star and occurs at values of the kinetic to 
potential energy ratio T/W ~ 0.27 [78]. For smaller 
rotation rates two secular instabilities with I = 2 arise, 
each corresponding to a sign of m = ±2. For incom- 
pressible fluids at constant density the T/W values for 
the onset of secular instabilities coincide. One insta- 
bility is driven by the viscosity, the other instability 
is driven by the gravitational radiation. For realistic 
stars the T/W values for the onset of these instabili- 
ties do not coincide; relativity and other factors shift 
the viscosity-driven instability to higher values of T/W. 
At the same time the gravitational radiation instability 
is shifted to lower values of T/W. The gravitational 
radiation instability arises for the modes which are ret- 
rograde in the co-rotating frame, while progradc in the 
(distant) laboratory frame. The underlying mechanism 
is the well established Chandrasckhar-Friedman-Schutz 
(CFS) mechanism [79, 80]. The bulk and shear viscosi- 
ties can prevent the development of the CFS instabil- 
ities, except in a certain window in the rotation and 
temperature plane. 

In the following we shall concentrate on axial modes 
of Newtonian stars, the so-called r-modes, which are 
known to undergo a CFS-type instability. Our main 
goal will be to assess the role of strong magnetic fields 
and bulk viscosity on the stability of these objects. We 
shall adopt the formalism of Refs. [22, 77, 81] for our 
study of the damping of the r-modes by bulk viscos- 
ity. For the sake of simplicity we shall describe both 
fluid mechanics and gravity in the Newtonian approxi- 
mation. 

The equations that describe the dynamical evolution 
of the star are 



d tP + v ■ Guv) = o, 



(79a) 



<9 f v + v- Vv= -V(/i-$) = -VU, (79b) 



V 2 $ = -A*kG P , 
where h is defined by the integral 



h(P) 



dP' 



(79c) 



(80) 



The quantity p is the mass density of the fluid which 
is assumed to satisfy a barotropic equation of state, 
p = p{P). $ is the gravitational potential and G is 
the gravitational constant. The potential U is used to 
determine the velocity field v. 

The oscillation modes of a uniformly rotating star 
can be completely described in terms of two perturba- 
tion potentials SU = U — Uo and 5$ = $ — $ , where 
Uq and $o are the potentials that correspond to the 
equilibrium configuration of the star. We assume that 
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the time and azimuthal angular dependence of any per- 
turbed quantity is described by oc e twi+Mre¥ ' , where m 
is an integer and Q is the frequency of the mode in lab- 
oratory frame. Let SI denote the rotation frequency of 
the star and w denote the frequency of the perturbed 
quantity measured in the co-rotating frame [which cor- 
responds the u> in Eqs. (48) and (49) because we will 
work in the co-rotating frame]. For small SI there is a 
simple relation between Sl,ui, and u> [77, 81], 



u = u) + mSl. 



(81) 



By linearizing the Euler equation (79b) around the 
equilibrium configuration, the velocity perturbation 
Sv a is determined by [77, 81] 



Sv a = iQ ab \7 b 5U. 



(82) 



The tensor Q ab is a function of u> and the rotation fre- 
quency SI of the star, 



Q ab = 



1 



(w + mSl) 2 - ASl 2 



(& + mSl)S 



ah 



An 2 



U) + mSl 



-z a z h 



2iV a v b 



,(83) 



where z is a unit vector pointing along the rotation 
axis of the equilibrium star, which we assume to be 
parallel to the magnetic field, i.e., = b{ in Cartesian 
coordinate system. Here Vo = rSl sin Otp is the fluid 
velocity of the equilibrium star. 

Having the linearized Euler equation, one proceeds to 
the linearization of the mass continuity equation (79a) 
and the equation for the gravitational potential (79c); 
one finds 



following discussion. The zeroth-order contribution to 
the r-mode is generated by the following form of the 
potential SUq, 



SU = 



m+l 



P™ +1 (cos0)e 



iu)t-\-im(p 



(to- 1)(to + 2) 



SI, 



(86) 



(87) 



where a is an arbitrary dimensionless constant and 
P[ n {x) are the associated Legendre polynomials. It has 
been shown that the most unstable mode is the one 
with to = 2 [22, 82], therefore we shall consider only 
this case in the following discussion. Substituting SUq 
into Eq. (82) one obtains the first-order perturbed ve- 
locity, 



Sv = a RSI 



-) Y* m (0,v>)e-', (88) 



a0r(m + 1) 3 (2to + l)!/m and Yf m (0,(p) 



where a' 

is the magnetic-type spherical harmonic function, 



r x VI/ 



I J! I 



(89) 



It is straightforward to check that the first-order per- 
turbed velocity satisfies 



d5vo z _ 
dz 

V-(5v = 0, 



(90) 



V a {pQ ab V b 5U) = -{Cj + mSl)(SU + 5<l>)dp/dh, 

V 2 <5$ = -4ttG(5U + 5$)dp/dh. (84) 

These equations, together with the appropriate bound- 
ary conditions at the surface of the star for SU and at 
infinity for 5$, determine the potentials 5U and 5$. 

For slowly rotating stars, Eq. (84) can be solved order 
by order in SI, 



SU = R 2 Sl 2 
(5$ = R 2 Sf 



SU - 
5$o 



su 2 



n 2 



irGp 
SI 2 
irGp 



0(S1 4 ) 

o(fi 4 ; 



(85) 



where R is the radius of nonrotating star. Since we need 
only the perturbed velocity, we will focus on SU in the 



therefore it does not contribute to the dissipation due 
to the bulk viscosities (± and <j|- In order to see how 
the bulk viscosities damp the r-mode instability one 
must consider next-to- first order, i.e., the third-order 
perturbed velocity which is generated by the poten- 
tial SU 2 - One can not determine analytically SUi from 
Eqs. (84)-(85), but the angular structure of 8U2 can be 
well represented by the following spherical harmonics 
expansion [81], 

SU 2 = «/ 1 (r)P^ +1 (cos0)e^ t+ ™^ 

+ah(r)P% +3 ( C os e)e iM+im *. (91) 

The functions fi(r) and f 2 (r) have been determined 
numerically in Ref. [81]. A useful approximation is pro- 
vided by the following simple expressions 



/i(r) = 
Mr) = 



-0.1294(4 
\R 



-0.0092 



(5)' 



0.0136 - 
\R 



0.1985 (I) -0.0388 (I 



0.0273 



R 



0.0024 (I) , 



(92) 
(93) 



which excellently fit the numerical result. We will use Eqs. (92) and (93) in the following numerical calcula- 
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tion. 

The energy of r-modes comes both from the ve- 
locity perturbation and the perturbation of the grav- 
itational potential. For slowly rotating stars, the 
main contribution comes from the velocity perturba- 
tion [22, 23, 82, 83]. Then, the energy of the r-mode 
measured in the co-rotating frame is 



Assuming spherical symmetry, we have 



(94) 



~ 2 



pr 2m+2 dr. (95) 



This energy will be dissipated both by gravitational 
radiation and by the thermodynamic transport in the 
fluid [22, 23], 



dE I dE 
~dt \~dt 



dE\ 

G \ / T 



(96) 



The dissipation rate due to gravitational radiation is 
given by [22, 23, 84] 



dE 
~dt 



where 



n)J2^ 21 [\SD lm \ 2 + \sj t 



l>2 



N, 



4irG{l + l){l + 2) 
Z(Z-l)[(2Z + l)!r 



(97) 



(98) 



SDi m and 5 Ji m are the mass and current multipole mo- 
ments of the perturbation, 



5D, 



I ill 



SJlr 



I 



r^v + tyvJ-Y^r. (99) 



Taking into account Eq. (87) one obtains 
2(m-l)(m + 2) o2 



Q(Q + mQ) = — - 



(m + 1) 



< 0, (100) 



which implies that the total sign of (dE/dt)c is pos- 
itive: gravitational radiation always increases the en- 
ergy of the r-modes. 

In order to compare the relative strengths of different 
dissipative processes, it is convenient to introduce the 
dissipativc timescalcs defined by 



2E 



(dE/dt)i 



(101) 



where the index i labels the dissipative process. 

The lowest-order contribution to {dE/dt)c comes 
from the current multipole moment SJu. For the most 
important case I = m = 2, this leads to the following 



timcscale (derived for a simple polytropic equation of 
state P oc p 2 ) [81, 82] 



1 



n 



i 



tq 3.26 \irGp 



2 \ 3 



(102) 



The bulk viscosities £j_ an d C|| dissipate the energy 
of the r-mode according to 




<Tv, 



(103) 



Accordingly, the time scales t^ ± and are given by 



2E 



T C-L,C|| 



(dE/dt) CxA 



(104) 



Currently, the shear viscosities j]± — 775 of strange 
quark matter are not known. In order to determine 
the damping of the r-modc by shear viscosity, we take 
as a crude estimate the value of 770 in the absence of a 
magnetic field [61] 



no 



r, = 5.5 x 10-V 5 / 3 /^ 473 x T" 5 / 3 , (105) 



where a s is the coupling constant of strong interaction. 
We will choose the value a s = 0.1 and apply Eq. (105) 
to highly degenerate 3-flavor quark matter with equal 
chemical potentials of all flavors, (p, u ~ p,d — Ms)- The 
contribution to the energy dissipation rate E due to 
shear viscosity r\ now becomes 



rilwij -%0/3|'cPr. (106) 



Assuming a uniform mass density star, the time scale 
t v can be simply expressed as [84] 




1 7?7 

~j ~ Tb 2 ' 



(107) 



The total time scale r(f2, T) is given by the following 
sum 



1 _ 1 1 1 1 



(108) 



which characterizes how fast the r-mode decays. Most 
importantly, if the sign of r is negative the amplitude of 
the r-mode will not decay, rather it will increase with 
time. Thus, it is important to determine the critical 
angular velocity f2 c for the onset of instability 



0. 



(109) 



At a given temperature, stars with f2 > il c will be 
unstable due to gravitational radiation. 

Figure. 8 shows the critical angular velocity fi c of a 
strange quark star with mass M = 1.4M Q and radius 
R = 10 km as a function of the magnetic field B near 
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FIG. 8: (Color online) The critical angular velocity f2 c of a 
strange quark star as a function of magnetic field B at tem- 
perature T = 0.001 MeV. The red dashed lines correspond 
to the unstable, while the blue solid lines correspond to the 
stable regime. 



10 18 G. The temperature is fixed as T = 0.001 MeV 
and other parameters are taken according to Eq. (52). 
In obtaining Fig. 8, we have taken into account the 
thcrmodynamical and hydrodynamical stability condi- 
tions. The blue-solid curves correspond to the thcrmo- 
dynamically and hydrodynamically stable region, while 
the red-dashed curves correspond to unstable regions. 
The critical angular velocity is strongly oscillating with 
increasing B. This behavior is due to the oscillating 
nature of the bulk viscosities 0i and (± as shown in 
Fig. 6. Thus, this macroscopic behavior originates from 
a purely quantum mechanical effect, namely the Lan- 
dau quantization of the energy levels of quarks. As 
discussed for the bulk viscosities G and (± shown in 
Fig. 7, averaging is needed to obtain physically rele- 
vant quantities. In Fig. 9 we show the averaged criti- 
cal angular velocity at various temperatures. The solid 
black curves are obtained by averaging over a short 
period Alog(B/G) = 0.05, whereas the short-dashed 
red curves correspond to averaging over a long period 
Alog(B/G) = 0.5. It is seen that after short-period av- 
eraging, the critical angular velocity (solid black curves) 
shows regular oscillation, the amplitude of which is 
growing as the £?-ncld increases. The critical angular 
velocity tt c displays a sharp drop for fields B < 10 18 5 G 
(short-dashed red curves), which is the consequence of 
the sharp drop of G and £j_ shown in Fig. 7. Thus, we 
conclude that for extremely large magnetic fields, the 
critical angular velocity at which the r-mode instability 
sets in could be significantly lower than in the absence 
of magnetic field. 

Figure. 10 shows the window of the r-modc instability 
in the 51 — log 10 T plane for a strange quark star of mass 
M = IAMq and radius R = 10 km. The regions above 
the respective curves correspond to the parameter space 
where the r-mode oscillations are unstable, i.e., a star in 
this region will rapidly spin down by emission of gravi- 
tational waves. The dashed green curve corresponds to 
vanishing bulk viscosities G| — (± =0. The solid black 
curve represents the (in) stability window of an un- 
magnetizcd strange quark star. The curves with sym- 
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FIG. 9: (Color online) The averaged critical angular veloc- 
ity fi c of a strange quark star as a function of magnetic 
field B at various temperatures T — 0.001, 0.1, and 10 
MeV. The black solid curve corresponds to averaging over 
a short period (Alog(B/G) = 0.05); while the red short- 
dashed curve corresponds to averaging over a long period 
(Alog(B/G) = 0.5). 



bols show the typical instability window for magnetic 
fields around 10 17 G (the red curve marked by triangles) 
and 10 18 8 G (the blue curve, marked by circles). The 
symbolled curves are obtained by using the bulk vis- 
cosities averaged over the period Alog 10 (-B/G) = 0.5. 
For low temperatures, T < 0.3 keV, the r-mode insta- 
bility is suppressed mainly by the shear viscosity; at 
these low temperatures the bulk viscosities are an in- 
significant source of damping, independent of how large 
the magnetic field is. However, for larger temperatures 
the bulk viscosities dominate the damping of r-mode 
oscillations. For magnetic fields below B <~ 10 17 G, 
the critical rotation frequency is almost independent of 
the S-ficld. The r-modc instability window increases as 
the magnetic field grows. For fields B > 10 18 G it is a 
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very sensitive function of the field, as a consequence of 
the rapid variation of the bulk viscosities with the field. 
Asymptotically, the instability window can become sig- 
nificantly larger than the window at zero magnetic field 
(see also Fig. 9). For completeness, Fig. 10 also shows 
the observed distribution of Low Mass X-ray Binaries 
(LMXBs) by the shadowed box, which corresponds to 
the typical temperatures (2 x 10 7 — 3 x 10 s K) and ro- 
tation frequencies (300-700 Hz) of the majority of ob- 
served LMXBs [85]. It is seen that even in the case of 
extremely large magnetic fields, our instability window 
is consistent with the current LMXB data. 



_._ b~10 188 G/ 


i 1 i i i 1 


\ B ~ 10 17 g/ 


— 1 












\ \ 
• V 


| LMXBs | 





U.Ui , 1 1 1 , 1 1 1 1 1 , I I I 

-5 -4 -3-2-10 1 2 
Log, (T/MeV) 

FIG. 10: (Color online) The r-mode instability window for a 
strange quark star. The star is stable below the respective 
curves. The green dashed curve corresponds to vanishing 
bulk viscosities Oj — (± — 0. The black solid curve rep- 
resents the window of unmagnetized strange quark matter. 
The curves with symbols show the typical behavior of the 
instability window when the magnetic fields are around 10 17 
G (red curve) and 10 18 ' 8 G (blue curve). The shadowed box 
represents typical temperatures (2 x 10 7 — 3 x 10 8 K) and ro- 
tation frequencies (300-700 Hz) of the majority of observed 
LMXBs [85]. 



VII. SUMMARY 



the fluid hydrodynamically unstable. 

For a number of reasons (density dependence of pa- 
rameters along the star profile, formation of domains, 
intrinsic multicomponent nature of the magnetic held) 
the dependence of the transport coefficients on the mag- 
netic fields are needed at different resolutions, i.e., they 
require some suitable averaging over a range of mag- 
netic field. We have provided such averages over an 
increasingly larger scale. We find that if the averaging 
period is small the bulk viscosities show regular oscilla- 
tions, the amplitudes of which increase with magnetic 
field (see the black solid curves in Fig. 7). These os- 
cillations are smoothed out if we further increase the 
averaging scale. At this larger scale the most interest- 
ing feature is the rapid drop in the bulk viscosity of the 
matter due to the confinement of quarks to the lowest 
Landau level; this occurs for magnetic fields in excess 
of B > 10 18 5 G (see the short-dashed red curves in 
Fig. 7). 

As an application, we utilized our computed 
anisotropic bulk viscosities to study the problem of 
damping of r-modc oscillations in rotating Newtonian 
stars. We find that the instability window increases as 
the magnetic field is increased above the value B > 10 17 
G. By increasing the field one covers the entire range of 
parameter space which lies between the two extremes: 
the case when bulk viscosity vanishes (dashed green 
curve in Fig. 10, which corresponds to extremely large 
magnetic fields B > 10 19 G for which the bulk vis- 
cosity drops to zero), and the case when the mag- 
netic field is absent (solid black curve in Fig. 10). The 
found novel dependence of the r-mode instability win- 
dow on the magnetic field may help to distinguish quark 
stars from ordinary neutron stars with strong magnetic 
fields, since the latter are much more difficult to magne- 
tize. It would be interesting to sec whether the objects 
that lie in between these extremes, e.g., hybrid config- 
urations featuring quark cores and hadronic envelopes 
(see ref. [86] and references therein), may interpolate 
smoothly between the physics of ordinary and strange 
compact objects. 



In this paper we have studied anisotropic hydrody- 
namics of strongly magnetized matter in compact stars. 
We find that there are in general eight viscosity co- 
efficients, six of them arc identified as shear viscosi- 
ties, the other two, £j_ and Q| , are bulk viscosities [see 
Eq. (23)]. We applied our formalism to magnetized 
strange quark matter and gave explicit expressions for 
the bulk viscosities [Eq. (48) and Eq. (49)] due to the 
non-leptonic weak reactions (la) and (lb). Due to the 
Landau quantization of the energy levels of charged 
particles in a strong magnetic field, the magnetic field 
dependence of £|| and £j_ is very complicated and ex- 
hibits "quasi-periodic" de Haas - van Alphen type os- 
cillations (see Fig. 6). For a magnetic field B < 10 17 G 
the effect of the magnetic field on the transport coeffi- 
cients is small and the bulk viscosities can be well ap- 
proximated by their zero-field values. For large fields 
10 17 < B < 10 20 G the viscosities are substantially 
modified, £j_ may even become negative for some val- 
ues of the B-field. We showed that negative £j_ render 
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Appendix A: Linear independence of components 
in Eq. (21) 



Contracting Eq. (Al) with and b" we find the fol- 
lowing two additional conditions 



The purpose of this appendix is to show explicitly 
that the eight different decompositions in Eq. (21) are 
linearly independent. To this end, let us write down a 
general linear combination of the eight different decom- 
positions, 



+ a2 (h) + h a 8 (viii) = 0. 



(Al) 



If (i) — (viii) are linearly independent, the coefficients 
a\ — a g should all vanish for any values of the vectors 
u M and V 1 . Firstly, it is easy to see that the compo- 
nents (vi) and (vii) are independent of the other com- 
ponents, because they have odd parity under reflection 
¥ L —> — b^, whereas the other six components have even 
parity under this transformation. Besides that it is ob- 
vious that (vi) and (vii) arc independent of each other. 
Therefore, we only need to treat the linear equation 



ffli(i) + • • • + 05 (v) + ag(viii) = 0. 



(A2) 



By contracting the indices \x and u, we obtain the fol- 
lowing three conditions 



3ai + 2a 2 - a 3 + 2a 8 = 0, 
3a 3 — a 4 + 4a 5 + 2a 8 = 0, 
3a! + 2a 2 — 12a3 + a 4 — 4a 5 = 0. 



(A3) 



ai + a 3 = 0, 
2a 2 — ct3 + ci4 — 4a5 = 0. 



(A4) 



Contracting the indices v and a in Eq. (61), we obtain 
one further condition 



<2l + 4d 2 — 2d3 + 04 — 605 = 0. 



(A5) 



The non-trivial solution of the set of Eqs. (A2)-(A5) is 



ai = (23 = 0, 
a 2 = a 4 /2 = a 5 = -ag. 

Then, we have the following condition, 



(A6) 



The only possible solution is a 2 = 0, which thus proves 
the independence of (i) — (viii). 
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